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Abstract We summarize briefly the main results obtained within the proposed 
EDABI method combining Exact Diagonalization of (parametrized) 
many-particle Hamiltonian with Ab Initio self-adjustment of the single- 
particle wave function in the correlated state of interacting electrons. 
The properties of nanoscopic chains and rings are discussed as a func- 
tion of their interatomic distance R and compared with those obtained 
by Bethe ansatz for infinite Hubbard chain. The concepts of renormal- 
ized orbitals, distribution function in momentum space, and of Hubbard 
splitting as applied to nanoscopic systems are emphasized. 



1. Introduction 

Recent development in computing techniques, as well as of analytical 
methods, has lead to a successful determination of electronic properties 
of semiconductors and metals based on LDA [1], LDA+U [2] and related 
[3] approaches. Even strongly correlated systems, such as V2O3 (under- 
going the Mott transition) and high-temperature superconductors, have 
been treated in that manner [4]. However, the discussion of the metal- 
insulator transition of the Mott-Hubbard type is not as yet possible in 
a systematic manner, particularly for low-dimensional systems. These 
difficulties are caused by the circumstance that the electron-electron in- 
teraction is comparable, if not stronger than the single-particle energy. 
In effect, the procedure starting from the single-particle picture (band 
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structure) and including subsequently the interaction via a local poten- 
tial might not be appropriate then. In this situation, one resorts to 
parametrized models of correlated electrons, where the single-particle 
and the interaction-induced aspects of the electronic states are treated 
on equal footing (in the Fock space though) [5]. The single particle wave 
functions are contained in the formal expressions for model parameters. 
We propose to combine the two efforts in an exact manner, at least for 
model systems. 

Our method of approach to the electronic states grew out of the 
following question: Can one complete the procedure starting from a 
parametrized model by determining the single-particle wave functions 
in the resultant correlated state a posteriori? In other words, we de- 
termine first the energy of interacting particles in terms of the micro- 
scopic parameters rigorously and only then optimize this energy with 
respect to the wave functions contained in those parameters by deriving 
the self- adjusted wave equation for this state. Physically, the last step 
amounts to allowing the single-particle wave functions to relax in the 
correlated state. This method has been overviewed in number of papers 
[6, 7, 8, 9, 10], so we present here examples of its application to low- 
dimensional and nanoscopic systems. We start with the analysis of the 
infinite Hubbard chain and then compare the results with those for finite 
chains. We also discuss briefly small hydrogenic ring of N = 6 atoms. 
Also, throughout the paper we are using adjustable Wannier composed 
of atomic or Gaussian functions, which are determined explicitly from 
the minimization of the system ground state energy as a function of 
interatomic distance. The paper describes various ground-state charac- 
teristics of simple monoatomic chains and rings. 

2. Exact diagonalization combined with 

wave-function determination: formal aspects 

As our method contains both many-particle and single-particle (wave- 
function) aspects, both treated in a rigorous manner, it may be useful 
to summarize the basic principle behind it. First, we start with the 
standard expression of the many-particle Hamiltonian in the Fock space 
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where Hi and H2 are the Hamiltonian for one and one pair of particles, 
and 



is the field operator, {wi(r)} is the single-particle basis of wave- functions 
(complete, but otherwise arbitrary), ad is the annihilation operator 
of the particle in the single-particle state \ia > represented by Wi(r) and 
the spin quantum number a = ±1. The only approximation we make 
in our whole analysis is that instead of taking the summation over a 
complete set {i} of single-particle states (and transition between them), 
we limit ourselves to a finite subset of M states. This means that we 
are solving a model many-body system rather than the complete problem 
(Hubbard or extended Hubbard models are classic examples representing 
one-orbital-per-atom) . 

The essential step in our analysis follows from taking the finite single- 
particle basis, which amounts to limiting the occupation- number repre- 
sentation space to a space of finite dimension. To minimize the error in 
estimating the ground-state energy of many-particle system we calculate 
first all configurations in the limited Fock subspace and then optimize the 
orbitals in the interacting (correlated) ground state. In this manner, the 
second quantization takes care of counting various many single-particle 
microconfigurations enforced by the interaction between them, whereas 
the wave- function optimization adjusts each of them to the milieu of 
all others. In brief, second-quantization aspect addresses the question 
how they are distributed among the single-particles state and the first- 
quantization optimization tells us how their states look like once they 
are there. 

One should also address the problem of many-body vs. single-particle 
wave function. In the wave mechanics of the interacting system only 
the iV-particle wave function \P(ri, . . . r^v) has a sense. However, in the 
second quantization the single-particle wave function appears explicitly 
in the expression for the field operator, the remaining part is the evalua- 
tion of various microconfigurations, with proper weights characterized by 
their energy (no entropy appears as they form a single coherent state). In 
other words, the microconfiguration counting in the occupation-number 
representation replaces the determination of TV-particle Hilbert space. 
But then, we are faced with the single-particle states determination, on 
which the counting is performed. Explicitly, the TV-particle state |$o > 
in the Fock space can be defined as 



$(r) = $>i(r) 




(2) 
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where |0) is the vacuum state. The iV-particle wave function is then 
determined from 

*o(ri,...,r N ) = ^=<0|*(n)...*(r N )|* >- (4) 

Expanding |<E>o > in the basis involving M single-particle states, i.e. 

-, M 

\^) = ~m E C n ..,J n ...a) N n, (5) 

V ' 31,— Jat=1 

we obtain the iV-particle wave function in the form 

j M M 

*o(ri,...r n ) = — E (°l«^---«n4--- a Ll°) 

' ii,...,iiv=l ji,—,3N=l 

Cj 1 ...j N w h (r 1 )...w iN (r N ). (6) 

The many-body coefficients Cj 1--- j N will be calculated from either the 
direct diagonalization or the Lanczos algorithm, whereas the wave func- 
tions {wi(r)} will be determined from the renormalized (self- adjusted) 
wave equation, which is set up in the following manner. First, we substi- 
tute (2) into (1) and obtain the formal expression for the ground state 
energy 

E G = (H) = J2 tij(4 a a ja ) + ^ Yl VijkMoA^i^k^), (7) 

ija ijkla\oi 

where the microscopic parameters 

t l0 = J cPrwtWHWwjir), 

and 

V ijk i = J d^r 1 d 2 d z r2W*{r 1 )w*{Y2)V{Y 1 -r2)w k {r 1 )wi{r 2 ), 

contain the single-particle wave functions and the averages are (a\ a aj a ) = 
(^olaj^Ojo-l^o), etc. Second, in the situation when we work with defi- 
nite number of particles (or else, if the chemical potential can be re- 
garded as constant), then we can determine {wi(r)} by setting Euler 
equations for each of them, i.e. by minimizing the functional F = 
Ec{wi(r), Vitfj(r) — J2i^if d 3 rw*(rwi(r. Such a procedure leads to the 
equation 

5E G ^ SE G 

V x/ _ = AiWi(r) (8) 



Sw*(r) <5(Vu£(r)) 



5 



Single-particle 
Schrodinger eq. 



£ ,iiijWj(r) = ejWi(r) 



Single-particle 
(or trial) basis 

Field operators 

Diagonalization 
in the Fock space 



Mr)} 



^(r),^ 



H = \*o)E G (%\ + . 



Ground-state energy E G = (*o| H |*o) 



Single-particle \ 
basis optimization J 



r en (r),(r en ) t (r) 



Renormalized N-particle 
wave function 

Figure 1. Flowchart of the EDABI method described in this paper. The top line is 
missing when the adjustable Gaussian basis is used. 



A direct solution of this equation is very difficult to achieve. Therefore, 
we define the starting wave functions wl°\r) = J2jLi Pij^j{ r i a )> where 
flij are the mixing coefficients and <3?j(r; a) are the atomic wave functions 
of the size a -1 . In result, the renormalized wave equation reduces to the 
minimization of (7) with respect to a (there is only one size a -1 when 
we take orbitals of the same type for each atomic site in the system). 
The whole EDABI procedure is schematically summarized in Fig. 1. 

3. Electron states for the Hubbard chain and a 
comparison with nanochains 

We implement first the EDABI method to the case of the linear chain 
composed of N atoms, which obeys periodic boundary conditions. In 
the simplest situation we have one valence electron per atom, as is in 
the case of monoatomic chains composed of Na, K or Cs. However, for 
the sake of simplicity, we consider here the chain composed of hydrogen 
atoms; this means only that we compose the Wannier function of the 
conduction band from ls-like atomic states of an adjustable size (there 
is no principal problem in considering ns-like states, with n > 1). 

As mentioned earlier, we start from the many-body model of inter- 
acting electrons. For this purpose, we consider an extended Hubbard 
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model, as represented by the Hamiltonian 



+U Y n i1 n U + Y K ij n i n 3 + 



(9) 



i i<j i<j 



The first term represents the atomic energy (e a = (wi\Hi\wi), where 
Hi is the Hamiltonian for a single particle in the system and Wi = 
Wi(r) is the Wannier state centered on site i in that system). The 
second term is the so-called hopping term with the hopping integral 
t = (wi\Hi\wi±i)(we use the tight-binding approximation and disregard 
more distant hoppings). The next two terms represent respectively the 
intra- and inter- atomic parts of the Coulomb interaction between elec- 
trons, with Kij = (wiWj\Vi2\wiWj) , U = Ka, and where Vi2(r — r') is 
the classical Coulomb interaction for pair of electrons. The last term 
expresses the classical repulsion between the ions located at sites i and 
j. Also, the single particle operator Hi contains both kinetic energy 
and the Coulomb attractive interaction of ions (it is sufficient to take 
6 10 ionic Coulomb wells surrounding the electron located on site i to 
reproduce to a good accuracy the effective values od e a and t). 

For a detailed analysis it is convenient to rewrite (9) in the equivalent 
form, which for the case of one electron per atom reads 



H = e e J f N e +tJ2 (4a a i+i<r + h.c^+U^2n i ^n ii +^2K ij Sn i 6n j , (10) 



where Srii = 1 — nj, N e = J2i a n «°- ^ s the number of electrons (here equal 
to the number of atoms) , and 



is the effective atomic energy, which includes the compensating repulsive 
interactions to provide correctly the atomic limit. We disregard the last 
term in (10), since we would like to relate the nanochain results to those 
for the Hubbard chain for N — > oo [11]. Under these circumstances, 
Hamiltonian (10), we consider here explicitly, represents only the Hub- 
bard model with the energy part e e J^ providing a proper neutral-atom 
limit for R — > oo. Eq. (10) then can be diagonalized exactly with the 
help of Bethe ansatz [12]. Explicitly, the expression for the ground state 
energy in terms of microscopic parameters e e J$ , t, and U takes the form 



i<j 




(11) 




(12) 
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where J n (x) is the Bessel function. One should note that the energy 
expression is not an additive function of terms ~ t and ~ U, as the 
solution is of nonperturbative nature. 

As we have stressed earlier, the Lieb-Wu solution (12) does not repre- 
sent the final step of the analysis as it still contains parameters, which in 
turn are expressed through the one-particle (Wannier) functions {wi(r)}. 
Therefore, we minimize the energy functional E = E{uii(r), X7wi(r)} 
with respect to {wi(r)}. Such a procedure leads to the Euler varia- 
tional principle for renormalized or self- adjusted wave functions in the 
correlated state [5]. Only after solving that equation and calculating 
explicitly the parameter values for given R we obtain the energy of the 
correlated ground state as a function of the lattice parameter. This last 
step completes the theoretical analysis of the Hubbard model. 

Technically, we compose the Wannier functions of atomic ls-like func- 
tions of variable size, with respect to which we minimize E. Explicitly, 
in the spirit of tight-binding approximation we can write that 

^(r) = /3^(r) - 7 [V>m(r) + &-i(r)], (13) 

where (3 and 7 are the mixing coefficients determined from the orthonor- 
mality condition (uii\wj) = Sij and 

^(r) = (a 3 /vr) 1/2 exp(-a|r - R,|), (14) 

with a being the adjustable parameter. Substituting (13) to the expres- 
sions for e a , t, and U, and then those expressions to (12) we obtain the 
physical ground state energy as the minimal energy with respect to a 
for given lattice constant R = \Ri — Ri±i\. 

In Fig. 2 we plot the optimized ground state energy with respect 
to a as a function of interatomic distance. We have also added there 
the corresponding values obtained using the Gutzwiller ansatz (GA)[13] 
and the Gutzwiller wave-function (GWF)[14] approximations. Obvi- 
ously, the approximate solutions [13, 14] represent the upper estimates 
for the system energy. In the inset we have plotted the interaction-to- 
bandwidth ratio U /W; one sees that the electrons are strongly correlated 
(i.e. U/W > 1) for realistic values of R (qq ~ 0.53 A is the Bohr ra- 
dius). Obviously, the chain composed of hydrogen atoms is not stable, 
as the energy Eq (per atom) is above — lRy. However, this is not an 
issue here, since we are considering only a model situation (such a chain 
could be stabilized on a substrate, but then we should have to include 
the trapping potential, in addition to the periodic potential composed 
of the protonic Coulomb wells). It would be interesting to repeat this 
analysis for 2s and 3s orbitals representing the conduction band of a 
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Figure 2. Optimized ground state energy as a function of interatomic distance 
- solid line, Gutzwiller ansatz(GA) solution - dotted line, Gutzwiller wave-function 
(GWF) - dashed line. 



quantum wire composed of Li and Na respectively; this does not present 
a major obstacle. 

In Fig. 3 we draw the Wannier function centered on site "0" for the 
hydrogen chain: the solid line represents the self-adjusted Wannier func- 
tion in the tight-binding approximation, whereas the dashed curve is the 
usual Wannier function calculated in TBA. For comparison, the dotted 
line is the usual ls-type atomic function. All the curves were drawn 
along the chain direction. Although the differences in the first two cases 
do not seem crucial, the values of microscopic parameters: e e J$ , t, U 
and K differ remarkably. The determined values of those parameters 
versus R/ao are provided in Table 1 (the values of Eq and a are also 
listed there). One notices numerically, that U/W > 1, where W = 4|i|, 
for R>2a . 

The detailed electronic properties of the chain have been discussed 
separately [15]. Probably, the most interesting result of relevance to this 
workshop is the conclusion that the values of the energy per atom (and 
of microscopic parameters as well) are almost the same when obtained 
either from the solution for N — > oo (discussed above) and from the 
numerical solution for, say, N = 10 atoms [16]. This statement is il- 
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Figure 3. The shape of the optimized Wannier function - solid line, non-optimized 
Wannier function - dashed line, Is atomic wave function - dotted line (for comparison). 



R/aa 


aminflO 




t 


U 


K 


E a /N 


1.5 


1.806 


0.9103 


-1.0405 


2.399 


1.695 


0.0665 


2.0 


1.491 


-0.1901 


-0.5339 


1.985 


1.172 


-0.5179 


2.5 


1.303 


-0.6242 


-0.3076 


1.722 


0.889 


-0.7627 


3.0 


1.189 


-0.8180 


-0.1904 


1.553 


0.713 


-0.8800 


3.5 


1.116 


-0.9104 


-0.1230 


1.440 


0.596 


-0.9391 


4.0 


1.069 


-0.9559 


-0.0815 


1.365 


0.513 


-0.9693 


4.5 


1.039 


-0.9784 


-0.0546 


1.317 


0.451 


-0.9848 


5.0 


1.022 


-0.9896 


-0.0370 


1.288 


0.403 


-0.9926 


6.0 


1.013 


-0.9977 


-0.0165 


1.269 


0.334 


-0.9982 


7.0 


1.001 


-0.9995 


-0.0072 


1.252 


0.286 


-0.9996 


8.0 


1.001 


-0.9999 


-0.0031 


1.251 


0.250 


-0.9999 


10.0 


1.000 


-1.0000 


0.0003 


1.250 


0.200 


-1.0000 



Table 1. Optimized inverse orbital size, microscopic parameters and the ground-state 
energy for N = 10 atoms calculated in Slater-type basis, as a function of interatomic 
distance. Intersite Coulomb repulsion K\ is included on the mean-field level in e%** , 
Hubbard U term is treated exactly. Single-particle potential contains six Coulomb 
wells. 
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E G 


R/ao 


N = oo 


N = oo 


N = 10 


N = oo 


N = oo 


N = 10 




ls/3G 


3G 


3G 


ls/3G 


3G 


3G 


1.5 


1.2985 


1.3062 


1.3094 


-0.5788 


-0.5527 


-0.5684 


2.0 


1.1753 


1.2038 


1.2047 


-0.8246 


-0.8104 


-0.8154 


2.5 


1.0924 


1.1203 


1.1203 


-0.9152 


-0.9123 


-0.9139 


3.0 


1.0485 


1.0683 


1.0672 


-0.9540 


-0.9560 


-0.9567 


4.0 


1.0212 


1.0212 


1.0203 


-0.9832 


-0.9840 


-0.9841 


5.0 


1.0109 


1.0062 


1.0054 


-0.9939 


-0.9900 


-0.9901 


6.0 


1.0055 


1.0020 


1.0027 


-0.9981 


-0.9914 


-0.9914 


7.0 


1.0021 


1.0005 


1.0000 


-0.9995 


-0.9917 


-0.9917 


8.0 


1.0005 


1.0001 


1.0000 


-0.9999 


-0.9917 


-0.9917 


10.0 


1.0002 


1.0001 


1.0000 


-1.0000 


-0.9917 


-0.9917 



Table 2. Renormalized values of the inverse size (a) of the atomic wave-function 
(columns 2 -f- 4) and the corresponding values of the ground state energy (columns 
5-5-7), both versus the lattice parameter R. 



lustrated in Table 2, where the optimal inverse size a m i n of the atomic 
wave functions composing iOj(r) and that of ground-state energy have 
been listed as a function the interatomic distance. The three columns 
in each category represent respectively the following results: (i) when 
single-site and two-site parameters have been calculated for ls-like func- 
tions and 3- and 4-site terms (in the atomic basis) have been calculated 
in the contracted STO-3G basis; (ii) and (hi) represent all calculations 
in the Gaussian basis. The results are pretty close, independently of 
the trial atomic basis selected to represent Wi(r). However, there is one 
restriction, namely the boundary conditions selected in the N — > oo and 
N = 10 cases must coincide ( here they are selected as periodic b.c). 
The importance of the results displayed in Table 2 relies on the fact 
that in this manner we can apply analytic results obtained for the infi- 
nite chain to the finite chains (nanowires) if only realistic single particle 
wave functions are taken into account. Obviously, in Table 2 the compar- 
ison was made for the case of one electron per atom only (the Lieb-Wu 
solution provides then the insulating state). We plan calculating the 
properties of Li and Na nanowires starting from this prescription. Par- 
enthetically, since in the infinite chains we have charge-spin separation 
and other non-Fermi liquid effects [17]; they should appear also in some 
form in non-half filled nanoscopic chains, which should also be regarded 
as strongly correlated systems from the start. This last question will be 
taken up again in Sec. 5. 
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4. Nanoscopic H N rings 

As a second example we consider a ring composed of N = 6 hydrogen 
atoms arranged planarly (the stable H4 clusters arranged spatially have 
been considered elsewhere [16, 5]). In this situation, the periodic bound- 
ary conditions (PBC) are the physical condition for the system geometry. 
In Fig. 4 we plot the profile of the wave function located around the 
exemplary atom in the hexagon Hq. This density contains renormalized 
orbitals and the calculations of the ground states involve (in principle) 
( 6 ) = 924 6-particle states in the occupation number representation 
spanned on M = 12 states and containing 6 Wannier functions of ad- 
justable size. The space profiles are useful for the determination of the 
density function profile n(r) = (^(r)^(r)), where ^(r) is the field op- 
erator spanned on those 6 Wannier states. In effect, we have that 



n(r) = \wi{r)\ 2 (n ia ) + Y^w*{v) Wj {v){a\ a a ja ), (15) 

icr ijcr 



where the primed summation means that i 7^ j. The second part rep- 
resents the part which does not appear in the Hartree-Fock (single de- 
terminant) approximation for the many particle wave function. In Fig. 
5 we display the electron density profiles (normalized to unity) for the 
N = 6 atoms arranged in the hexagon, as is also in the translationally 
invariant along the ring density profile n(r). 

An interesting feature of the spectrum of electronic states arises 
when the distance between the atoms increases. Namely, the spectrum 
decomposes into well defined Hubbard subbands, as shown in Fig. 6 
(more appropriately, they represent manifolds corresponding to the sub- 
bands when N ^ 00). The lowest manifold (I) corresponds to the con- 
figuration with approximately singly occupied orbitals (highest occupied 
Wannier orbitals) whereas the manifolds II-IV correspond respectively 
to the states with one to three double occupancies. This division into 
the well separated manifolds for larger R is even better seen for the 
clusters of N = 4 and 5 atoms [8]. One should mention that the states 
considered here represent the excited states calculated with the help of 
Lanczos procedure [5] , repeated many times until the configuration with 
the minimal energy and the optimal single-particle wave function size 
are reached simultaneously. 
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Figure 4- Spatial wave-function profiles for a selected site of He cluster of hexagonal 
shape. Note the negative value on the neighboring site to the central atom. 
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Figure 5. Exact density profiles for electrons in a hexagonal ring of atoms. 
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Figure 6. Decomposition of the system energies into Hubbard subbands for H§ 
cluster, plotted as a function of interatomic distance. The horizontal lines represent 
the atomic limit values for the levels S, with I = 0, 1,2, and 3. 



5. Further features of the results: collective 
properties of nanochains 

5.1 One electron per atom case: localization 
threshold 

In the previous Sections we illustrated the applications of the ED- 
ABI method, in which the interaction among particles is dealt with first. 
This is because, in most cases, the interaction parameters (coupling con- 
stants) represent the largest energy scale in the system. The first of the 
examples (the Hubbard chain) represents the situation, for which an an- 
alytic expression for the ground state energy exists [12], whereas the 
case of Hn rings must be treated numerically all the way through [18]. 
In applications of this method to the extended three-dimensional sys- 
tems one will have to resort to the approximate treatments of the model 
Hamiltonian in the Fock space. This last problem poses a real challenge 
for the future. In the remaining part of this brief review we concentrate 
on the collective properties of the nanochain. 

We have concentrated first on the basic quantum-mechanical fea- 
tures of the system such as the ground-state energy or the renormal- 
ized single-particle wave function in the milieu of other particles. In 
Fig. 7 we plot the exact ground state energy of a chain of N = 6 -f- 10 
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atoms and compare it with that obtained in the Hartree-Fock approxi- 
mation for the Slater antiferromagnetic state. The starting Hamiltonian 
is of the form (10). As one can see the Hartree-Fock energy repre- 
sents an upper estimate, as it should be. Additionally, the curve M 
represents the "metallic" approximation, for which the correlation func- 
tion < Sriidrij > has been taken for the ID electron gas on the lat- 
tice. On the contrary, INS represents the energy of the Heisenberg- 
Mott state in the mean-field approximation. The state of the sys- 
tem crosses over from the Slater metallic-type state to the localized- 
spin-type of state. This is seen explicitly when we calculate the evo- 
lution of the spin-spin correlation function with the increasing inter- 
atomic distance, as displayed in Fig. 8. Well defined oscillations of 

< Sj • Sj > are seen for even (N = 12) number of atoms, which become 
more pronounced with the increasing N (the frustration effects appear 
for odd N). What is much more important, the autocorrelation part 

< Si-Si >=< Sf >= (3/4)(l— 2 < n^nn >) evolves from the value close 
to the free-electron value < Sf >= (3/4) (1 - 2 < X >) = 3/8 
to the atomic-limit value < Sf >= (l/2)(l/2 + 1) = 3/4. This evolution 
provides a direct evidence of the crossover from delocalized to the local- 
ized regime. 

Other properties such as the electrical conductivity [16] and the sta- 
tistical distribution in momentum space (rika) have also been addressed 
[9]. Here the question emerges whether the quantum nano-liquid of elec- 
trons in a nanochain resembles at all the Landau-Fermi liquid or if it is 
rather represented by the Tomonaga-Luttinger scaling laws [17]. The an- 
swer is not yet settled, as within our method we can deal only with up to 
N=16 hydrogen atoms assembled into a linear chain. However, one can 
make some definite statements for the particular cases. Namely, for the 
half-filled case (one electron per atom) the modified Fermi distribution is 
a good representation of the nk a for smaller R values; with the increasing 
atom spacing it is smeared out above critical spacing R = R c « 3.4ao 
[10], which characterizes a crossover from the case with extended states 
to the state of localized electrons on atoms, as detailed below. In Fig. 9 
we exhibit this evolution on the example of the distribution function (i.e. 
electron occupation in the momentum space); a clear universality is ob- 
served for N = 6 -=-14 atoms. The adjustable Gaussian (STO-3G) single- 
particle basis has been used in the analysis. The existence of the Fermi 
ridge quasi-discontinuity for a small R is very suggestive in this case and 
is positioned near the Fermi wave- vector = ±n/(2R), corresponding 
to that in Landau-Fermi liquid, with TV" — > oo. For N = 6, 10, 14 periodic 
boundary conditions (PBC) provide the minimal ground state energy, 
whereas for TV = 8 and 12 the anti-periodic boundary conditions (ABC) 
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a/ao 


DU 


#1*2 


Dlo 


D* 8 


D% 




a 




1.5 


0.9225 


0.9420 


0.9563 


0.9727 


0.9822 


0.8008 





019 


1.6 


0.8879 


0.9162 


0.9378 


0.9612 


0.9754 


0.7175 


o 


029 


1.7 


0.8419 


0.8817 


0.9130 


0.9459 


0.9667 


0.6148 





043 


1.8 


0.7826 


0.8365 


0.8805 


0.9256 


0.9552 


0.4967 





064 


1.9 


0.7095 


0.7794 


0.8389 


0.8992 


0.9406 


0.3728 





092 


2.0 


0.6245 


0.7105 


0.7875 


0.8660 


0.9222 


0.2567 





129 


2.1 


0.5315 


0.6310 


0.7265 


0.8254 


0.8996 


0.1606 





172 


2.2 


0.4338 


0.5431 


0.6549 


0.7755 


0.8714 


0.0899 





228 


2.3 


0.3403 


0.4523 


0.5766 


0.7179 


0.8379 


0.0455 





287 


2.4 


0.2554 


0.3631 


0.4937 


0.6524 


0.7982 


0.0207 





352 


2.5 


0.1839 


0.2812 


0.4109 


0.5813 


0.7526 


0.0087 





420 


2.6 


0.1269 


0.2096 


0.3315 


0.5065 


0.7009 


0.0033 





489 


2.7 


0.0840 


0.1508 


0.2595 


0.4312 


0.6441 


0.0012 





566 


2.8 


0.0536 


0.1049 


0.1972 


0.3586 


0.5836 


0.0004 





641 


2.9 


0.0315 


0.0706 


0.1456 


0.2914 


0.5208 


0.0001 





920 


3.0 


0.0196 


0.0461 


0.1047 


0.2314 


0.4575 


0.0000 







Table 3. Normalized Drude weight D* N , the extrapolated value D^, and its relative 
error for ID half-filled system with long-range Coulomb interactions. 




2 3 4 5 6 7 

INTERATOMIC DISTANCE, R (a.u.) 



Figure 7. Ground state energy per atom vs. R for the linear chain with N = 
6-t-IO atoms with periodic boundary conditions. The STO-3G Gaussian basis for 
representation of atomic orbitals forming the Wannier function has been used. The 
inset provides a universal behavior of the inverse size a of the orbitals. For details 
see main text. 



16 




Figure 8. Spin-spin correlation function (S; • Sj) vs. the distance \i — j\ between 
the atomic sites, for different lattice constant R and for N = 12 atoms. Due to 
periodic boundary conditions only the distance up to \i — j\ = 6 is relevant. The 
continuous line is guide to the eye. A quasi-antiferromagnetic arrangement is clearly 
seen, particularly for larger R. 
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Figure 9. Evolution of the statistical momentum distribution rik a with the increas- 
ing interatomic distance R. For the nanochain of N = 6, 10, and 14 the periodic 
boundary conditions provide the energy minimum, whereas for N = 8 and 12 an- 
tiperiodic boundary conditions are appropriate. For details see main text. 
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INTERATOMIC DISTANCE, a (a.u.) 



Figure 10. Luttinger-liquid scaling for a half-filled ID chain of N = 6 -f- 14 atoms 
with long-range Coulomb interactions: (a) momentum distribution for electrons in the 
linear and (b) log-log scale, continuous lines represent the fitted singular expansion 
in powers of \n(n/\kF — k\a) (see main text for details); (c) Tomonaga-Luttingcr 
model exponent 9 vs. lattice parameter a (specified in ao) and (in the inset) the 
corresponding residual sum of squares. The solid lines in Figures (a) and (b) represent 
the TLM fitting of Eq. (16). 
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lead to the lower energy. However, one clearly sees the absence of the 
points at the Fermi momentum kp°. This is because, for example, the 
Fermi points for N = 14 are located at kpRjix = ±3/7, whereas they 
are located at kpRj-n = ±5/12 for N = 12, close to the values ±1/2 in 
both cases. One would have to apply a renormalization group approach 
[19] for the states close to Uf (i.e. perform the analysis for larger number 
N ~ 10 2 atoms) to determine the precise evolution of the distribution 
function, this time with the system size N. Nevertheless, the results for 
N < 14 atoms represent those for a true nanoscopic system. 

Before addressing the question of localization directly, we would like 
to address the question whether the computed distribution displayed in 
Fig. 9 can be fitted into the Tomonaga-Luttinger mode, with the log- 
arithmic scaling corrections included [17]. The statistical distribution 
near the Fermi point can be represented by 



with z = 7r/\k — Here 9 is a non-universal (interaction-dependent) 
exponent (it yields the nonexistence of fermionic quasiparticles, since its 
residue vanishes as ~ \k — kp\ 6 with k — ► kp. The corresponding 
electron-momentum distribution is depicted in Fig. 10a in the linear, 
and in Fig. 10b in the log-log scale. The R dependence of the exponent 
9 is shown in Fig. 10c and crosses the value 9 = 1 for R = R c m 2.6ao 
corresponding to the localization threshold [17]. This threshold is about 
30% smaller than the corresponding value (R c « 3.4ao) when the almost- 
localized Fermi-liquid view was taken [10]. One should also note that 
the Luttinger scaling does not reproduce well the occupancies farther 
away either way from the Fermi point. Thus the results concerning n^ a 
do not provide a definite answer as to the exact nature of the nanoliquid 
composed of N ~ 10 electrons, although it is absolutely amazing that 
they have such a nice and simple scaling properties. 

To address the question of the electron localization directly, we have 
calculated the optical conductivity a(uj), which can be written in the 
form a(uj) = D5(uj) + a reg (uj), where the regular part is 



In \nka- — np\ = — 9 In 2 + 6 In In z + 0(1/ In z), 



(16) 



Vreg{u>)— £ 



5{u - (E n - Eq)), 



(17) 



n^0 



E n — Eq 



whereas the Drude weight (the charge stiffness) D is given by 



£ = -<*o|T|* >-^£ 



< *n|jp|*0 > 



(18) 



n^0 



E n — Eq 
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with T being the hopping term as in (9) and j p the current operator de- 
fined as j p = itJ2ij C r( a ja a i^ ~ h.c). Here |* n > is the system eigenstate 
corresponding to the eigenvalue E n . For a finite system of N atoms D is 
always nonzero due to nonzero tunnelling rate through a potential bar- 
rier of finite width. Because of that, the finite-size scaling with 1/N — > 
must be performed on D. We use the following parabolic extrapolation 

lnD* N = a + b(l/N) +c(l/N) 2 ,} (19) 

where D* N = —{N/ir)D/ < ^ol^l^o > denotes the normalized Drude 
weight for the system of N sites. We observe that < D* < 1 and hence 
can be regarded as an order parameter for the transition to the localized 
(atomic) states. In Table 3 we plot the weights D* N , the extrapolated 
values -D£o, and its relative error for ID half- filled system with long- 
range Coulomb interactions included. What is very important, the value 
of drops by two orders of magnitude when R changes by a factor 
of two (between 1.5ao and 3ao). Note also that is within its error 
for R ~ R c ~ 2.3ao, close to the value obtained from the Tomonaga- 
Luttinger scaling for n^, as one would expect. The result for D* N is 
probably telling us how far we can go quantitatively when discussing the 
localization in nanowires. These results will be detailed in a separate 
publication. 

5.2 Quarter-filled case 

For the quarter-filled case (i.e. when every second atom contributes 
a valence electron) the distribution is more smeared out, as shown 
in Fig. 11. The ground state for R > 3ao is then well represented by 
the charge-density-wave state, since then the density-density correlation 
function ((nj — n)(rij — n)) exhibits the oscillatory behavior, as shown 
in Fig. 12. The onset of the charge-density wave state in that case is 
invariably due to the long-range Coulomb interaction ~ Kij, which is 
reduced in that state. The charge-density wave order parameter defined 
as 

Qcdw = ^Y.i- 1 )" 1 < (ni-n)(nj-n) > (20) 

■m 

reaches its maximal value 1/4 for R > 8ao- Let us stress again, the form 
of the statistical distribution (and its feasibility) for iV ~ 10 atoms is a 
very interesting question by itself, since this is the regime of nanoscience. 
Our results show that even in that regime one should be able to see the 
signatures of the phase transition to the spin- or charge- density wave 
states. 




Figure 11. Momentum distribution n^a for electrons on a chain of N = 8-f- 16 atoms 
in the quarter-filled band case (N e = N/2). Lines are drawn as a guide to the eye 
only. Values of the lattice parameter R are specified in units of do- PBC and ABC 
denote periodic and antiperiodic boundary conditions, respectively. The dashed line 
marks occupation n ka = 1/2. 




Figure 12. Charge-density distribution for the quarter-filled (N e = N/2) nanochain 
as viewed by density-density fluctuation correlation function (An; An.,) vs. relative 
distance \i — j\. 



21 



6. Conclusions 

The EDABI method (Exact Diagonalization combined withAB Initio 
orbital readjustment) provides the exact ground state energy of the model 
systems considered (Hubbard chain, nanoscopic chains and rings) as a 
function of interatomic distance. It also provides reliably other ground- 
state dynamical characteristics for nanoscopic systems: spin and charge 
correlation functions, the spectral density and the density of states, as 
well as the system conductivity. Not all the characteristics have been 
presented in this overview [5, 10, 16, 17]. Furthermore, the method 
can also be extended to nonzero temperatures. Finally, it should be 
underlined again that our method of approach is particularly suited for 
strongly correlated systems, where the interaction and the single-particle 
parts should be treated on the same footing. The exact numerical results 
in a model situation can also serve as a test for approximate analytic 
treatments. 

The analysis with inclusion of long-range Coulomb interactions of the 
distribution function rika in either Fermi-like or Tomonaga-Luttinger 
categories suggests the existence of the crossover transition to the local- 
ized state with the increasing interatomic distance R. This state is either 
the spin- or the charge- ordered for the number of electrons N e = N and 
N/2, respectively. In the small R limit the state can be rendered as 
quasi-metallic in the sense that the quasimomentum hk can be regarded 
as a good quantum number, even though the level structure is discrete. 
The difference between the short-chain and infinite-chain situation is due 
to the circumstance that in order to form extended states in the present 
case the electrons tunnel through a barrier of finite width (the length 
L = NR). This is one of the reasons for quasi-metallicity. The other is 
the presence of the long-range Coulomb interaction [20]. 
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